Searching for the continuum spectrum photons correlated to the 130 GeV 7-ray line 
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Indications for a 7-ray line(s) signal towards the Galactic center at an energy of about 130 GeV 
have been recently presented. While dark matter annihilations are a viable candidate for this signal, 
it is generally expected that such a flux would be correlated to a 7-ray component with continuum 
energy spectrum due to dark matter pair annihilating into other Standard Model particles. We use 
the 7-ray data from the inner 10° x 10° window to derive limits for a variety of DM annihilation 
final states. Extending the window of observation, we discuss bounds on the morphological shape of 
a dark matter signal associated to the line, applying both standard templates for the dark matter 
profile, such as an Einasto or a NFW profile, and introducing a new more general parametrization. 



I. INTRODUCTION 

Dark Matter (DM) accounts for approximately 85% 
of the matter density of the Universe. Among the many 
possible scenarios on its nature, Weakly Interacting Mas- 
sive Particles (WIMPs) thermally produced in the early 
Universe are compelling candidates. WIMPs have a very 
rich phenomenology: they may be produced at colliders 
as the LHC, they may be detected by their direct inter- 
actions with baryonic matter or indirectly singling out a 
component in cosmic rays and 7-rays due to their pair 
annihilation into Standard Model (SM) particles. 

Among the possible signals of DM annihilation/decay, 
the search for a monochromatic 7-ray flux in multiple 
directions is one of the indirect detection methods with 
cleanest signature. As a prime target for such a signal, 
there is the galactic center (GC) direction, given that the 
annihilation/decay rate of DM in the halo peaks in that 
region. 

Recently [l[ has suggested the detection of a line at 
129.8 ± 2.4+ 7 GeV with a 3.3 a significance in a wide 
window towards the GC. Additionally has revealed 
an excess of 7-rays concentrated around the GC that, if 
modeled as a single line, is at 127.0±2.0 GeV at 5.0 a sig- 
nificance. A pair of lines with energies of 110.8 ± 4.4 and 
128.8 ±2.7 GeV can also explain the excess of 7-rays with 
5.4 a significance. Apart from the morphological differ- 
ences these results agree on the energy of the line. Also 
both signals arc in agreement with the constraints on line 
searches from the Fermi Collaboration Q . The morphol- 
ogy of the line signal from both [l[ and especially from Q 
leads us to concentrate on DM annihilation rather than 
decay. 

Most WIMP models that can give a linc(s) signal do so 
with a production rate which is loop-suppressed. Typi- 
cally the monochromatic 7-ray yield comes together with 
a 7-ray yield with continuum spectrum due to the annihi- 



lation, at tree-level, into other two-body SM final states, 
which in turn hadronizc and/or decay into the stable 
species, p, p e^, i^s and 7-rays. Given that we have a 
preferred direction (inner few degrees of the GC) and a 
preferred mass range (correlated to the ~130 GeV line), 
for the DM explanation of the line(s), we can place spe- 
cific limits on the continuous 7-ray component from DM 
annihilation (see for example [3, l3] ) ■ Alternative expla- 
nations for the line signal have also been suggested in 
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In section |TT] we describe our method of calculating 
limits on the continuous component for a variety of SM 
annihilation products. Additionally, using the 7-ray data 
in the energy range of the line(s) and assuming it is a DM 
signal we can study the morphology of events at that 
energy in the sky by extending our observation window. 
We show these results in section IIIII In section IIVI we 
implement our limits for a specific DM scenario and give 
our conclusions in section [V] 



II. EXTRACTING LIMITS ON DM 
ANNIHILATION FROM THE INNER 10° x 10° 

In Q the angular extension of the 7-ray line(s) signal 
at 110-130 GeV is described by a gaussian with FWHM 
of 3° —4°. If that signal comes from DM annihilation, the 
same region must be used to study the room for a 7-ray 
yield with continuum spectrum due to DM annihilation 
into other final states piling up at energies below 130 
GeV. For that purpose we use the region of the inner 
10° x 10° box (| b \< 5°, I I \< 5°). 

In this analysis we have calculated spectra using 3 
years of Fermi LAT 7-ray data, taken between August 
2008 and August 2011. Using the FermiTools PASS7 
(v9r23pl) "ULTRACLEAN" class of data that ensures 
minimal CR contamination [29j . We bin all the 7-ray 
events with energies of 200 MeV and up to 500 GeV 
in 30 logarithmically spaced energy bins [30j . We sepa- 
rately calculate the exposures (and fluxes using HEALPix 
0) for front and back-converted events, before summing 
their contribution to the total flux within the windows of 
interest. 
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To account for the contribution from known point and 
extended sources within the windows of interest we have 
used the 2 yr published catalogue (see [l(| and references 
therein). Since we care to extract conservative limits on 
DM annihilation from 7-rays with energy below 130 GeV 
towards the inner 10° x 10° we ignore the contribution 
of the Fermi Bubbles / Fermi haze [lH, (T^j; given 
that there are significant uncertainties on the exact mor- 
phology of these structures at lower latitudes [l3| . 

The 7-rays in the window of the inner 10° x 10° orig- 
inate from a combination of sources. There are 29 de- 
tected point sources centered in that window (To| 2 close 
by extended sources that contribute minimally [3, EH , 
as well as the diffuse 7-rays from inelastic collisions 
of CRs with the interstellar medium (ISM) gas, from 
bremsstrahlung radiation off CR electrons and from up- 
scattering of low energy photons of the interstellar radia- 
tion field (ISRF) from high energy CR e (inverse Comp- 
ton scattering). Additionally many unknown dim point 
sources are expected to be located within that window. 
Finally the possible 7-ray contribution from DM annihi- 
lations in the halo is expected to peak towards the galac- 
tic center. 

These 7-rays can be the direct product of DM annihi- 
lations as in the case of the monochromatic yield from 
the 27, Z-f or hn final states, possibly matching the lines 
detected by 0, . Also virtual internal bremsstrahlung 
(VIB) and final state radiation (FSR) in DM annihila- 
tion can give a very hard spectrum that can be confused 
as a line over an otherwise featureless power law spec- 
trum 

0,01. The decay of mesons (predominantly 7r°s), 
produced in the decay or hadronization processes of the 
products of DM annihilation, can also lead to a significant 
contribution to the gamma-ray spectrum. This compo- 
nent though while typically harder than the background 
7-ray spectra, is significantly softer than the VIB/FSR 
and can not be confused as a DM line in 7-rays. These 
contributions probe directly the DM annihilation pro- 
file and will be referred to, as prompt 7-rays. Addition- 
ally, inverse Compton and bremsstrahlung 7-rays from 
the lcptonic final products (e ± ) of DM annihilation will 
also contribute in that window since both the ISRF en- 
ergy density and the ISM gas density peak in the inner 
part of the Galaxy. 

As we suggested in the introduction, we will concen- 
trate only on the DM annihilation case since the line(s) 
morphology is so confined that it favors profiles of cuspy 
annihilating DM halos (see also our discussion in section 

EH). 

To compute the diffuse 7-ray background we use the 
DRAGON package [UGH [H with a new ISM gas model 
[20I ] that ensures goo d ag reement with 7-ray data in that 
window and overall [21|. We ignore in this work the 
contribution form the dark gas whose uncertainties are 
though significant in the inner 5° in latitude [22|, [23| but 
that would only add to the diffuse 7-ray background re- 
sulting in less room for DM annihilation originated 7- 
rays. 



We study five individual modes/channels of DM an- 
nihilation: XX — > W + W~, XX — > bb, XX — > t + t~, 
XX — ► and XX — ► e + e~. Typically, DM mod- 

els have sizable branching ratios into more than one of 
these channels. The exact limits in such models can be 
recovered by linearly combing the limits from the above 
channels. Annihilation to Z gauge bosons (top quarks) 
give very similar 7-ray spectra to those of W + W~ (b 
quarks) and thus the constraints to those channels can 
be taken to be the same (within ~ 10%) to those of the 
XX — ► W+W- (xx — > bb). 

Following we assume that a line at energy of 127 ±2 
GeV has been detected. The morphology of the excess 
is described by a bi-gaussian with FWHM of 4 degrees 
in both I and b. That line can come from XX — * 27 
or XX — ► -^7 or XX — ► ^7- in [3 a single line at 
129.8 ±2.4;!^ 3 GeV has been suggested. Additionally the 
case where there are 2 lines centered at 128.8 ± 2.7 and 
110.8 ± 4.4. GeV has been indicated by 0. In that case 
the lines come from either the combination of 27&Z7 
lines or from the Z7&/17 lines. 

We study both the case of a single line centered at 
127 GeV and the case of 2 lines centered at 129 and 111 
GeV. The choice of mass depends on the exact origin of 
the line(s). For a single line from XX — ^ 27 the mass 
range of 122 < m x < 132 GeV is studied. For a single 
line from XX — > %1 we study 137 < m x < 145 GeV 
and from XX — ► ^7 we study the 149 < m x < 157 
GeV mass range. For 2 lines originating from 27&Z7 
we study 127 < m x < 130 and for the case of Z7&/17 
lines 138 < m x < 143 GeV. The relevant ratio in the 
luminosity of the two lines is taken to be 0.7/1 for the 
111/129 GeV lines. We allow for 4% uncertainty in the 
determination of energy of the line(s) which leads to the 
ranges of masses referred above, which is about 2er of the 
declared uncertainties of and 0] [3l| . 

For every choice of annihilation channel to continuum 
7-rays, annihilation channel(s) to line(s) and DM mass, 
we first find the best fit values from the 7-ray data within 
10° x 10° for both the cross-section of the main annihi- 
lation channel (giving the continuum 7-rays) and sepa- 
rately for the annihilation cross-section to the line(s) (2 
d.o.f.). 

In Fig. [H we show a fit to the total 7-ray spectrum 
within our window of interest for the case of XX — > 
W + W~ , with m x = 130 and a single line coming from 
XX — ► 27. The fact that the best fit value for the 
cross-section is positive validates our claim of deriving 
conservative limits on DM annihilation, while the good 
agreement to the low (E 1 < 1 GeV) energies where the 
DM contributes at the few % level shows the good agree- 
ment of the physical model for the background to the low 
energy data. 

From the best fit value we then derive the 3a upper 
limits of the main annihilation channel keeping the an- 
nihilation to the line fixed to its best fit value. In Fig. [3] 
we show these 3<r upper limits for the five annihilation 
channels to W + W~, bb, t + t~, /i + ^ and e + e~. We 
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FIG. 1: Case of m x = 130 GeV DM particle annihilating to a 
W + W~ pair with a cross-section of 1.05 x 10 -25 ci 
to a 2y line with a cross-section of 1.25 x 10 -27 cr 
plot the | 6 |< 5°, | I \< 5°. 
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show limits for both the case of a single line Fig. [2] (top 
left) and for a double line Fig. [3] (top right). The exact 
choice of the origin of the linc(s) and its encrgy(ies) has 
a subdominant effect on the limits in all channels except 
in the case of e + e~. That happens since the DM con- 
tribution of the main channel to the 7-ray spectrum is 
at energies bellow 100 GeV where the lines do not con- 
tribute. The case of the e + e~ channel is an exception 
due to the very significant FSR component which peaks 
at m x . Thus the FSR component competes with the 
linc(s) in the fit, making its limits sensitive to the ex- 
act assumptions on the line(s). We also give in Fig. [2] 
(bottom panels) the best fit values for the line(s) for the 
relevant combinations of DM mass and channel. 

The ISRF photon and gas densities have been fixed 
based on our background model. The assumptions 
on these densities influence the inverse Compton and 
bremsstrahlung components respectively. One can derive 
even more conservative limits on the DM annihilation 
channels by considering only the prompt 7-ray contribu- 
tion. 

In Fig. [3] we give the 3 tr limits where only the prompt 
7-rays from DM are taken into account. For the W + W~ , 
bb and t + t~ channels, for which the prompt 7-rays are 
the dominant component, the limits become weaker only 
by ~ 10 — 20%. For the p + p~, e + e~ modes on the 
contrary, since hard CR electrons are injected, their in- 
verse Compton and bremsstrahlung components are sig- 
nificant. Thus if we ignore these diffuse components keep- 
ing only the prompt component, the 3er limits become 
weaker by a factor of 4-5 in both channels. 

The limits shown in Figs. [2] and [3] depend on the DM 
profile assumptions. We use here an Einasto DM profile: 



p(r) = pEin exp 



2 fr a 



1 



(1) 



the local DM density is equal to 0.4 GeV cm" 3 [U El- 
That results in a J- factor from that window of J/ AD. = 
1.21 x 10 24 GeV 2 cm -5 , where J factor is defined here as: 



J = 



An Jo 



p 2 DM ( s ,tydsdn, 



(2) 



with a = 0.22, R c = 15.7 kpc and pEin is set such that 



with s to be the distance along line of sight and Aft the 
angle of observation. 

A more cuspy DM profile would lead to stronger limits 
while a more cored (flat) in the inner kpes would lead to 
weaker limits. All the limits shown in Fig.[3]and the lim- 
its for W + W~ , bb and t + t~ in Fig. [2] will change inverse 
proportionally (exactly or approximately) to the value of 
the J-factors within that window, since the prompt com- 
ponent is dominant in these channels. The same applies 
for the best fit values to the line(s). Thus these limits 
can be used for other DM profile assumptions once one 
properly takes into account the different J-factor from 
that window. For the annihilation channels into p + pT 
and e + e~ the limits in Fig. [2] have a dependence on the 
DM profile that is more involved. 

Finally since our aim in this paper is not to study the 
line itself but the accompanying 7-ray fluxes for the DM 
case, we want to ensure that the exact line assumptions 
that we make do not influence our limits for the continu- 
ous component. The 3er limits presented in Figs|2]and[3] 
were derived with the cross-section to the line(s) to be 
the best fit value from the fit to the 7-ray data within 
I I |< 5°, I b |< 5°. Alternatively, we calculate the 3tr 
limits for the same channels using for the cross-section to 
the line(s) such a value that gives the luminosity stated 
for the 4° FWHM cusp of 0, that is (1.7±0.4) x 10 36 ph/s 
or (3.2 ±0.6) x 10 35 erg/s. The difference in the values of 
the cross-sections to the line(s) between the two methods 
is ~ 30% (at the same level with the stated uncertainty 
of i). 

In Table|I]we present our limits on the continuous com- 
ponents for these two alternative methods of evaluating 
the cross-section to the line(s) before deriving the limits. 
For the case where the cross-section value to the line(s) 
comes from the | / |< 5°, | b |< 5° region fit, (denoted 
as "free") and for the case where the cross-section comes 
form the luminosity stated by @] ■ We show all five chan- 
nels for three masses characteristic for the three DM mass 
ranges valid in the case of a single line at 127 GeV. 

The difference in the limits for all five channels and all 
masses between the two methods is at the ~ 1% level. 
The same results apply for the case of 2 lines (111 and 
129 GeV). Thus the exact luminosity assumptions for the 
line(s) can not influence our results on the continuous 
component. 

We also find that changing our window of observation 
from (| I |,| b |) < 5°, to (| I |,| b |) < 3°, 4° or 8° 
our limits for the continuous component (for the best fit 
value for the line(s)) change by up to 10% (20%), with 
the limits from (| I |, | b |) < 5° being the strongest (see 
also work of [2£ 
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FIG. 2: Top: 3 a limits for annihilation into "channel", from region of b |< 5°, | I |< 5°. We study five channels. 
XX — * W + W~: solid lines, XX — ► bb'- dotted lines, XX — > t + t~: dashed lines, XX — ^ / i+ A t_: dashed dotted lines and 
XX — ► e + e": long dashed lines. Left: assuming a single line from XX. — > 27 (blue), or XX — ^ %1 (red) or XX — ► ^7 (green). 
Right: assuming double lines from XX — ► 27 and XX — * -^7 (blue), from XX — ► ^7 an d XX — ► h"y (red). Bottom: best fit 
values for the ahhinilation into line(s). For the double line case the annihilation best fit value refers to the cross-section for 
the highest energy line; the 111/129 GeV luminosity ratio is taken to be 0.7/1 . We use the Einasto DM profile of eq. [T]which 
gives J I AQ = 1.21 x 10 24 GeV 2 cm" 5 (see text for more details). 



Chan. 


Line 


127 GeV (2 7 ) 


140 GeV (Z-y) 


150 GeV (hy) 


W + W~ 


Free 


34.2(40.8) 


35.1(42.6) 


36.6(44.1) 


W + W~ 


Fixed 


34.5(41.4) 


35.4(43.2) 


37.2(44.7) 


bb 


Free 


30.0(31.5) 


31.5(33.3) 


32.7(34.5) 


bb 


Fixed 


30.3(31.8) 


31.8(33.6) 


33.0(34.8) 


T + T- 


Free 


20.4(21.9) 


21.6(23.4) 


24.1(24.9) 


T+T- 


Fixed 


20.7(21.9) 


21.9(23.7) 


23.4(25.2) 


(j, + H~ 


Free 


39.0(155.7) 


39.9(169.8) 


42.0(185.4) 


(j, + H~ 


Fixed 


41.1(156.3) 


40.2(167.7) 


42.3(184.5) 


e + e~ 


Free 


18.3(91.8) 


13.5(100.8) 


18.9(111.0) 


e + e~ 


Fixed 


18.3(92.1) 


13.5(99.3) 


19.2(110.4) 



TABLE I: 3<r upper limits on DM annihilation (av) x BR 
to the continuous part (channel) (xlO -26 cm 3 s -1 ) using 
full(only prompt) DM 7 -ray spectra within | I |< 5°, | b < 5°. 
The line signal is taken to be either from its best fit value 
(free) or fixed using the luminosity of 0] (see text for more 
details). The J-factor/AJl from this window is 1.21 x 10 24 
GeV 2 cm" 5 . 



III. DARK MATTER ANNIHILATION SIGNAL 
PROFILE 

Ref. 0, suggests that the line(s) signal can be mor- 
phologically fitted by a 4° FWHM gaussian distribution 
or 3° FWHM when using just the events and avoid mak- 



ing diffuse maps or masking out any part of the GC. The 
author of [l| has suggested instead a wider region of best 
significance for the 130 GeV line. 

Using ULTRACLEAN data class we address as well the 
matter of DM profile morphology under the assumption 
that the line signal is of DM origin and that an associated 
continuous spectrum exists. 

Calculating the 7-ray spectral data within a wider re- 
gion of the sky we can derive limits on the allowed an- 
nihilation cross-section for a specific assumption on the 
DM halo profile or vice versa for the DM halo profile 
properties for specific assumptions on the annihilation 
cross-section. 

Motivated by the ~ 130 GeV energy of the 7-ray line, 
we consider a DM mass of m x = 130 GeV annihilating 
to W + W~ , with a cross-section to 2j for the line. 

We calculate the 7-ray spectra in the same energy bin- 
ning as for the 10° x 10° box described in section [TTJ We 
concentrate in the | b \< 25°, | / |< 25° region where the 
annihilation from the halo is dominant. We break that 
region in 20 smaller windows, 8 of which in the region 
I b |< 10°, I I \< 10° of 5° size in I and 10° size in b. For 
I b |> 10°, I I |> 10°, 12 windows symmetrically placed 
with respect to b = 0° , I = 0° each composed of 4 boxes 
5° x 5° size. These windows are also shown in Fig. 2J 

In Fig. |31 we assume the DM profile to be cither an 
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FIG. 3: Top: 3 a limits for annihilation into "channel" using only the prompt 7-rays from DM annihilation. We use the region 
of I b |< 5°, I / |< 5°. Left: assuming a single line. Right: assuming double lines. Lines and colors as in Fig. [2] Bottom: best 
fit values for the annihilation into line(s). We use the Einasto DM profile of eq. [1] 
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FIG. 4: Strength of limits on the annihilation cross-section to W W . Left: Einasto profile. Right: NFW profile. Values are 
normalized to the 3 a limit on the annihilation cross-section of \X — ► W + W~ from —5° < I < 0°, | b < 5° ((cru)^o <i<o° |b|<5° 
= 3.72 x 10 - 25 cm 3 s -1 ). Values smaller than 1 indicate stronger (lower) 3 a limit by the given (multiplication) factor. 



Einasto as described in eq. [TJ (Fig. HJ left panel) or an 
NFW profile (Fig.£Q right panel): 



P(r) = Pnfw 




(3) 



with R' c = 14.8 kpc, p N FW = 0.569 GcV cm- 3 . For both 
cases we fit the DM density to the locally measured value 
of 0.4 GeV cm -3 [24|. That results in a specific J-factor 



for each window. 

We calculate as in section ILTl the 3a limits on \X 
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W + W~ from each angular window. These limits include 
7-ray background contribution that as in the 10° x 10° 
box fits the lower energies. 

In Fig. S] we give the ratios of the 3er limit on each 
window devided by the 3a limit from the window of 
-5° < I < 0°, | b |< 5°: 

\ /window (a\ 

\ u u /-5 o <l<0°, \b\<5° 

Values of the ratio in eq. [4] smaller than 1 indicate 
stronger 3cr limits than that in the window of reference. 
We choose the window of -5° < I < 0°, | b |< 5° as 
reference since -as also thoroughly discussed in 0]- the 
majority of photons with energy ~ 130 GcV come from 
that part. We find 38 photons with energy 104.5 — 135.7 
GeV (the one energy bin in our analysis that includes 
the line(s)) to come from the -5° < I < 0°, | b |< 5° 
window, while only 21 photons in the same energy bin 
come from the symmetrical to the GC, 0° < / < 5°, 
\ b |< 5° window. 

In our simulations the DM halo profile is centered at 
I = 0°, 6 = 0°. We find that our 3 a limits are stronger 
from the 0° < I < 5°, | b |< 5° window than from the 
-5° < I < 0°, | b < 5° one. We remind that these limits 
come from the continuous component (see discussion in 
section [II]), thus there is an indication even from the con- 
tinuous DM component that there is an excess of 7-ray 
fluxes from the I < side. That is seen with both the 
Einasto and NFW profiles. 

We do not claim that such an excess is due to DM 
annihilations since it could equally correspond to an un- 
derestimation of the background in that energy range. 
Though that can be the case also for the line @, Q. Yet, 
if the line is of DM origin and is off-center as suggested 
by [1| , the strength of the derived limits may indicate an 
analogous effect in the continuum term. 

Once we move to | I |> 5° and | b |< 5° the strength 
of the 3cr limits drops. For the Einasto profile the lim- 
its become stronger at -10° < I < 5°, 5° <| b |< 10° 
and weaker in all other windows. For the more cen- 
trally peaked NFW profile the only region where slightly 
stronger limits are recovered is that of —5° < I < 0°, 
5° <| b |< 10°. 

Turning the perspective around, we discuss now what 
information can be extracted on the dark matter profile, 
having fixed the annihilation cross section to a reference 
value. Rather than introducing for the dark matter den- 
sity some functional form specified in terms of few param- 
eters, as most usually done in the literature, we consider 
here a much more general approach. Our generic spheri- 
cally symmetric profile is set via: i) Specifying its value 
Pi at the seven Galactocentric distances r.; = Rq sin(ai), 
being R the Sun Galactocentric distance and with on = 
5°,10 o ,15 o ,20 o ,25 o ,45 o &65°, namely at the spherical 
shell corresponding to the angular windows already in- 
troduced above, plus three higher latitude patches; ii) 
Fixing the dark matter density at the local Galactocen- 



tric distance r§ = Rq to p§ = 0.4GeVcm -3 and imple- 
menting a linear interpolation in a double logarithmic 
scale to retrieve the density profile between any two of 
these radii, i.e. allowing for an arbitrary power-law scal- 
ing between any two r^, with the only extra assumption 
of imposing that the profile is monotonically increasing 
for decreasing radius; Hi) Assuming that the profile fol- 
lows our reference Einasto model for r > Rq, a choice 
that has no impact on the analysis that follows. 

We refer again to the sample dark matter case intro- 
duced in Fig. [TJ with m x = 130 GeV and annihilating to 
W + W~ . While the ratio between the cross section into 
27 to the one into W + W~ is fixed by the best fit value 
(in the specific case we found it to be about 0.012), the 
absolute value of the cross section scales with the inverse 
of the line of sight integration factor J in the angular 
window I b |< 5°, | / |< 5°, in turn depending on the 
density profile within all the angular shells introduced 
above. After choosing a reference value for (trw), we wish 
to derive how large a contribution to J may come from 
each of the shells in our model, without violating the 
constraints coming from data in other angular windows. 
This gives an indication on how centrally concentrated 
the dark matter profile should be to provide a signal in 
the GC direction and, at the same time, to be consistent 
with data away from it. For this purpose we introduce the 
factors Ji which are analogous to the J-factor introduced 
in Eq. [2 except for imposing that the density profile is 
constant below the radius r^, namely p(r < n) = p\. 

The analysis is performed scanning the parameter 
space defined by the values pi (with i £ [1,7]). For each 
model wc compute line of sight integration factors cor- 
responding to the angular regions displayed in Fig. 21 as 
well as for the regions at 25° <| b |< 45°, 45° <| b |< 65° 
and 65° <| b |< 85°, and 0° | I |< 20° (for all three lat- 
itude intervals). For all regions we can compare against 
the 3cr upper bound on the flux due to a dark matter can- 
didate with mass 130 GeV and annihilating to W + W~ 
(conservatively including prompt emission only, while ig- 
noring the radiative emission from the associated lepton 
yields). As upper bound to the monochromatic signal 
we consider instead the mean flux integrated in the en- 
ergy bin [104.5, 135.7] GeV, obtained again using the UL- 
TRACLEAN data sample, and under the very conserva- 
tive hypothesis of zero background from diffuse emission. 
Among models passing constraints, we search for config- 
urations giving the maximum for the individual terms Ji 
(we also implement the additional limit Ji < J, given 
that any additional contribution to the line of sight inte- 
gral at radii r < r\ is always neglected in our setup). Al- 
though we are dealing with a very large parameter space, 
finding the upper bounds to J^, which we label J|" aa: , 
is not exceeding expensive since one can show that, for 
each radial shell, they mostly correspond to the models 
with largest changes in profile slope between neighboring 
shells. 

Ratios between J™ ax and J are shown in Fig. 03 where 
we display separately the J™ ax found when applying the 
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FIG. 5: Upper bounds on the partial Ji factor, normalized to 
the total GC line of sight integration factor J, for the para- 
metric dark matter density profile introduced in the text. The 
Ji factor are computed assuming a constant dark matter den- 
sity within the corresponding radial shell Vi. The J™ ax values 
displays are derived implementing separately the limits on the 
monochromatic flux (lower green bands) and the continuum 
spectrum (higher orange bands) derived from the other an- 
gular windows considered in the analysis. Limits are shown 
as narrow bands since they refer to three different values of 
the annihilation cross section, see the text for details. Also 
shown are the values for Ji/ J for our reference Einasto profile 
(blue dots) and for a cored Burkert profile (red squares) ; as it 
can be seen the Einasto profile is allowed, while the Burkert 
shape is excluded. 



limit on the monochromatic flux and when implementing 
that from the component with continuum spectrum; lim- 
its are shown as (very narrow) bands since they were 
derived for three different values for (o~v): the "ther- 
mal" value 3 x 1CP 26 cm 3 s^ 1 , the best fit value in case 
of our reference Einasto profile 1.05 x 10~ 25 cm 3 s -1 , 
and ten times the thermal value (this shows that depen- 
dence of our analysis on (o~v) is really very mild). For 
comparison, we plot also values of Ji/J for our refer- 
ence Einasto profile and for a Burkert profile, namely 
p oc l/(r + R c )/ (r 2 + i? 2 ), with local dark matter density 
0.4 GeV cm~ 3 and core radius R c = 10 kpc [13|. As one 
can see the Burkert profile is excluded from both line and 
continuum components, while the line limits are giving 
stronger evidence towards the need for a more centrally 
concentrated dark matter profile. This is most probably 
related to the fact that the limits arc derived in part from 
regions of the sky where the Fermi Bubbles/haze, has 
been claimed to be needed; we do not try to include such 
component in our background model and most probably 
this translates into an extra room (or a less severe con- 
straint) on the continuum emission from dark matter an- 
nihilations. On the other hand, the Fermi Bubbles/haze 
are expected to play a marginal role at high energy, hence 
the sharper constraint from the line cmissivity. 
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FIG. 6: Wino/Axion model of |2g|. m x 

{ov) xx — >z ~ l = 1.26 x 1(T 26 cmV 1 , {ovf° x = 3.2 x 10" 
cm 3 s~ 1 . 



IV. A SPECIFIC EXAMPLE 

The limits that we show in Fig. [2] and [3] can be used 
(linearly combined) for a wide class of models. As a spe- 
cific example we use our code to evaluate the constraints 
on a relevant model presented for the explanation of a 
130 GeV line 

In [28| the DM is composed by Winos and Axions 
at about equal amounts in DM mass density towards 
the GC. Following [28[ we take the DM mass density 
in Winos to be 49% of the total in the GC and in the 
entire Galaxy. Using the Einasto model of eq. Q] we take 
DM mass to be m x = 145 GeV and the cross-section 
to Z7 line to be 1.26 x 10 - 26 . The total annihilation 
cross-section of the Winos is 3.2 x 10 -24 cm 3 s _1 and is 
dominantly to W + W~ . Assuming a BR=0.96 for annihi- 
lation to W + W~ we derive that such a model is excluded 
as we show in Fig.[6l In fact such a cross-section is O(10) 
larger than the relevant 3tr limit for that mass and chan- 
nel shown in Fig. [5] 

We note that even ignoring the inverse Compton and 
bremsstrahlung components and all the background con- 
tribution, the prompt component which also includes the 
line signal overshoots the total j-ray spectrum between 10 
and 40 GeV. This result can not depend on just a differ- 
ent assumption for the DM profile or on a varying ratio 
in Wino to Axion mass density within the Galaxy since 
by changing any of these assumptions the 7-ray line will 
decrease/increase by the same amount. 



V. CONCLUSIONS 

Inspired by recent indications for a monochromatic 7- 
ray signal towards the GC, possibly connected to DM 
annihilations, we have derived limits on the continuous 
component that in general accompanies such signal. We 
use the region at | I |< 5°, | b |< 5° and compute 3a 
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upper limits for DM annihilation cross-sections into the 
W + W~, bb, t + t~ , fJ, + fi~ and e + e~ channels. 

We derive limits both from the total DM 7-ray emis- 
sion, including the prompt, the inverse Compton and the 
bremsstrahlung components as given in Fig. [2 and only 
from prompt DM 7-rays (given in Fig. [3]). These limits 
do not depend on the exact normalization of the linc(s) 
since they are dominated by the 7-ray data below 100 
GeV, where the lines do not contribute (see Table UJ). 
While our limits depend on the choice of the DM halo 
profile, they can be easily rescaled to a different halo 
configuration. This happens since, apart from the 
and e + e~ channels, the prompt 7-rays are the dominant 
component. 

We study the 7-ray data from other angular windows 
of the Galaxy and find that for cuspy DM profiles such 
as the NFW profile, the most stringent constraints come 
from our | I |< 5°, | b |< 5° window while for the 
Einasto profile slightly stronger limits can come from 
5° <| b |< 10° (see Fig. g]). We have also introduced 
a new general parametrization for the DM profile to dis- 
cuss how centrally concentrated the profile should be to 



give a flux compatible with the suggested GC line signal 
without violating bounds from other angular windows. 
We produced results for partial line-of-sight integration 
factors which are readily applicable to any dark matter 
profile, showing e.g. that a Burkert DM halo cannot be 
compatible (see Fig. 0). 

Our limits on specific DM annihilation channels can be 
linearly combined and readily applied to most DM mod- 
els in this mass range. We apply our limits to the model 
of [28| and conclude that it is excluded given that its 
prompt component exceeds the total 7-ray flux (Fig. [5]). 
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